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Abstract 

We investigate numerically 10 - dimensional Euclidean quantum 
gravity (with discretized Einstein - Hilbert action) in the framework of 
the dynamical triangulation approach. For the considered values of the 
gravitational coupling we observed two phases, the behavior of which 
is found to be similar to that of the crumpled and elongated phases of 
3, 4 and 5 dimensional models. Surprisingly, (for the observed lattice 
sizes) the natural state of the 10 D system (when the Einstein - Hilbert 
action is turned off) is found to resemble branched polymer while in 
the low dimensional systems the natural state belongs to the crumpled 
phase. 



1 Introduction 



The attempts to unify fundamental interactions have produced a lot of vari- 
ous models (see, for example, pQ and references therein). All of them are of 
rich mathematical structure and most of them are to some degree based on 
the Riemannian geometry. That's why we believe that the quantization of 
the latter is rather important. Even if it has nothing to do with the real grav- 
ity, it may play an important role in the further construction of the realistic 
unified theory. 

Recently it has been paid much attention to the dynamical triangulation 
(DT) approach to quantization of Riemannian geometry j2] in 2, 3 and 4 
dimensions 1 . For D = 2 the DT model has a well - defined continuum limit 
consistent with the predictions of the continuum theory j^j. At D = 3,4 the 
Euclidean DT models with the discretized Einstein action have two phases: 
the crumpled phase with infinite fractal dimension and the elongated one, 
which resembles branched polymer model with the fractal dimension close to 
2. It appears that the introducing of the causal structure to these models 
(correspondent to the transition from Euclidean to Lorentzian quantum grav- 
ity) or coupling them to matter changes their behavior and make them more 
realistic E] . Nevertheless, the pure Euclidean gravity is still of interest 
as an area of developing the investigation methods. Moreover, mechanisms 
observed within this model can be related to more realistic models. 

The mathematical structures related to the unification of fundamental 
forces may include the concept of high dimensional (D > 4) space. In par- 
ticular, all known superstring models become unambiguous only in the 10 
- dimensional space - time. Unfortunately, most of information about their 
structure comes from the perturbative methods and the (incomplete) inves- 
tigation of certain special excitations (such as D - branes) . It seems that the 
nonperturbative investigation of those models would become exhaustive only 
if the numerical lattice methods are used. The modern lattice theory has an 
experience of dealing with quantum gravity models in lower dimensions. In 
particular, the dynamical triangulation method was applied to the systems 
of D = 2,3,4 dimensions. Therefore we guess that it is reasonable to apply 
this method to the higher dimensional problems. The investigation of the 
pure Einstein Euclidean gravity can become the first step. The work related 
to this step is the content of the present paper. 

1 5 - dimensional Euclidean quantum gravity was considered in [3]. 
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Namely, we consider the ten - dimensional Euclidean dynamical triangu- 
lations of spherical topology. Their behavior was expected to be similar to 
that of 3,4 and 5 - dimensional systems. This supposition was confirmed 
partially. However, it turns out that there are a few features that are not 
present in the lower dimensional models. For example, the 10D model does 
have two phases that are similar to the crumpled and elongated phases of 3, 4 
and 5 - dimensional models. But in our case the phase transition (at least for 
the observed volumes V = 8000 and V = 32000) is at negative gravitational 
coupling constant. This means that at "physical" positive gravitational con- 
stant G the considered 10 - dimensional system cannot exist in the crumpled 
phase unlike 3, 4 and 5 - dimensional systems. Also this means that for the 
observed lattice sizes the natural state of the 10 D system (when the Einstein 
- Hilbert action is turned off) resembles branched polymer while in the low 
dimensional systems the natural state belongs to the crumpled phase. 

It is worth mentioning that the considered volumes are not consistent with 
the system of "physical" dimension D = 10. It follows from the observation 
that already at linear size 3 the analogous rectangular lOD lattice has the 
volume V = 3 10 = 59049. However, the observed effective (fractal) dimension 
of the triangulated space, say, at V = 32000 is ~ 4 in the crumpled phase 
and ~ 2 in the elongated one. So, the linear size of the system is expected 
to be of the order of ~ (32000) a ~ 10 and ~ (32000)^ ~ 100 respectively. 
This is confirmed by the direct measurement of the linear extent. 

2 The model. 

In this section we shall remind briefly the definition of the considered model 
and the description of the numerical algorithm. For the complete review 
of the dynamical triangulations see [2] and references therein. For the full 
description of the numerical algorithm for arbitrary dimension see [7j. 

In the dynamical triangulation approach the Riemannian manifold (of 
Euclidian signature) is approximated by the simplicial complex obtained by 
gluing the D - dimensional simplices. Each simplex has D + 1 vertices. All 
links are assumed to have the same length a. The metrics inside each simplex 
is supposed to be flat. Therefore the deviation from the flatness is concen- 
trated on the boundaries of the simplices. The scalar curvature R is zero 
everywhere except the bones (D — 2 - dimensional subsimplices of the trian- 
gulation simplices). The Einstein - Hilbert action can be expressed through 
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the number of bones and the number of simplices of the given triangulated 
manyfold: 



S = 



^ J R(x)^d D x = E (2* - 0(bone)cos-\±)) 



bones 



Vob_ 2 £>(£> + !) Ar m 

\-l 'bones , ^* simplices COS y ) ) , 



where O(bone) is the number of simplices sharing the given bone, Volj = 
a3 jW^ ls ^ e v °l ume of a i - dimensional simplex with the edges of length 
a, Agones is the total number of bones and Ar simplices is the total number of 
simplices. 

The metric of the triangulated manyfold is completely defined by the way 
the simplices are glued together. Therefore, the functional integral over Dg in 
this approach is changed by the summation over the different triangulations 
(we restrict ourselves with the spherical topology only): 

/^-Et^ (2) 

Where the sum is over the triangulations T that approximate different Rie- 
mannian manyfolds 2 and Ct is the symmetry factor of the triangulation itself 
(the order of its automorphism group). 

We consider the model, in which the fluctuations of the global invariant 
D - volume are suppressed. The partition function has the form: 



(3) 



Z V = J2 7^exp(-S(T)) = E 7 ^exp(/ tD _ 2 iV D _ 2 - k d N d - 7 (iV D - Vf 

where we denoted N D = N simp n ces , N D _ 2 = N honcs , and k d _ 2 = V °gg" 2 • 
Unfortunately, it is not possible to construct an algorithm that generates the 
sum over the triangulations of the same volume. So the constant 7 is kept 
finite. kd(V, is chosen in such a way that 

< N D >= E -£-exp(K D - 2 N D _ 2 - k d N d - j(N D - Vf)N D = V (4) 

T 



2 Two formally different triangulations may approximate the same Riemannian many- 
fold. Therefore the correspondent factorization should be implemented. 
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This provides that the volume fluctuates around the required value V. The 
fluctuations are of the order of SV ~ -^=. In order to approach the model 

with constant Nd we must keep ^ << 1. In practice we use 7 = 0.005 for 
V = 16000,32000. So ~ 10~ 3 . 

For the numerical investigations we used Metropolis algorithm in its form 
described in j7j. It is based on the following Markov chain. Each step of the 
chain is the proposal of a deformation T — > Tf of the given triangulation T iy 
which is accepted or rejected with the probability ViTi — > Tf) that satisfies 
the detailed balance condition 

eM-S{mV{T % - T f ) = eM-S(Tf))V(T f - T t ) (5) 

The definition of the proposed deformations is based on the following 
idea. Let us consider some closed D - dimensional simplicial manyfold of 
the topology of a D - dimensional sphere. Then, if a connected piece of our 
original triangulation is equal to a piece of this manyfold, we can replace it 
by the remaining part of the given manyfold. Thus we obtain the deformed 
triangulation with the same topology as the original one. Further, let us 
choose the boundary dso+i of a D + 1 - dimensional simplex sd+i as the 
mentioned above manyfold. There are D + 1 opportunities to distinguish a 
piece of dso+i correspondent to p - dimensional subsimplices of sd+i (p = 
0, D). The resulting deformation is called (p, D — p) move [7j. 

It has been shown that via such moves it is possible starting from any 
triangulation to reach a triangulation that is combinatorially equivalent 3 to 
an arbitrarily chosen other triangulation [S]. This property is called ergodic- 
ity Due to this property starting from the triangulation that approximates 
jH] the given Riemannian manyfold it is possible to reach a triangulation 
that approximates almost any other Riemannian manyfold. The exceptional 
cases, in which this is impossible, are commonly believed not to affect phys- 
ical results. 

In practise we choose randomly the type of the move (p G {0, D}), the 
simplex of a triangulation and its p - dimensional subsimplex. After that we 
check is there a vicinity of this p - dimensional simplex that is equivalent to 
the required piece of ds^+i- If so, the suggested move (and the correspondent 

3 Two triangulations are combinatorially equivalent if they have subdivisions that are 
equivalent to each other up to the relabeling of (sub)simplices. The subdivision T s of the 
given triangulation T is another triangulation such that the set of its vertices contains all 
vertices of T and any simplex s € T s either belongs to T or belongs to some simplex of T. 
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subsimplex) is called legal and we proceed with checking the possibility to 
perform the move 4 . If the move is allowed we accept or reject it with the 
following probability: 

P<T ' "* T,) = i+d+ ^y' i^TO-sm)) (6) 

In our calculations we start from the triangulation of minimal size that 
is dso+x- Then we allow it to grow randomly until the volume reaches the 
given value V. After that we proceed with the normal Metropolis process. 
During this process the coupling is self - tuned automatically in order 
to get its required value satisfying (JJJ). This is done via redefining it as 
k, — > k + 27(< Nd > —V) after each 10 sweeps 5 . 

In order to make a check of our results, we have made two independent 
program codes. The main program was written in C ++ using modern meth- 
ods of object - oriented programming. The algorithm is partially based on the 
one described in [7j and on the ideas suggested in [TU]. The second program 
was written in Fortran 77 and was used for checking the main program at 
small lattice sizes. The calculations were made within the parallel program- 
ming environment using the computation facilities of Joint Supercomputer 
Center (Moscow). 



3 Numerical results. 

We investigate the behavior of 10D model for V = 8000 and V = 32000. We 
considered the values of kd varying from —0.1 to 0.1, where the system is 
found to exist in two phases. The phase transition point is at kd ~ —0.03 

at V = 8000 and at k d 0.01 at V = 32000. 

The self - tuned value of k w appears to be independent of V. In accor- 
dance with large volume asymptotics obtained in jH] its dependence upon k, 8 
is linear with a good accuracy. The best fit is: 

kio = 15.57(1)k 8 + 0.42(1) (7) 

We investigate the following variables, which reflect the properties of the 
triangulated manyfold. 

4 The move is not geometrically allowed if the resulting new simplex already exists in 
the given triangulation. 

5 One sweep is V suggestions of legal moves. 
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1. The mean curvature carried by a simplex. The normalization is chosen 
in such a way that it is defined as 



R = 47rAW es 1 (8) 

D(D + lJJVrimpliceBOOS-Hi) 

The results for V = 8000 and V = 32000 coincide with a good accuracy 
They are represented in the figure ^ in contrast with the four dimensional 
case (at the observed values of « 8 ) R > (in the four dimensional case it 
becomes negative at the value of k^ 2 close to 0). The best fit to the curve 
of the figure [T] is: 

R = 0.218(1) + 0.654(2)k 8 - 1.97(3)«jj (9) 

2. The geodesic distance between two simplices is the length of the short- 
est path that connects them. The points of this path correspond to the 
simplices. The links of the path correspond to the pairs of neighboring sim- 
plices. We denote the geodesic distance between the simplices u and v by 
p(u,v). Let us fix a simplex s. Then the ball Br(s) of radius R G Z is 
consisted of the simplices u such that p(u, s) is less than or equal to R. The 
volume of the ball is defined as the number of simplices contained in it. 

One of the most informative characteristics of the triangulated manyfold 
is the average volume V(R) of the balls (of the constant radius R) as a 
function of this radius. In practise during each measurement we calculate 
Vr(s) for arbitrary chosen simplex s. Then we perform its averaging over the 
measurements (separately for each R). V(R) becomes constant at some value 
of R. This value is the averaged largest distance between two simplices of the 
manyfold, which is also called the diameter d. We found that for V = 8000 
at 1 < R < | the dependence of log V on log R is linear. For V = 32000 
the same takes place at 4 < R < |. The slopes give us the definition of the 
fractal (Hausdorff) dimension of the manyfold: log V = const + D\ogR. 

3. The mean fractal dimension T>(V) of the manyfold for the observed 
volumes V as a function of kd-2 is represented on the Fig. 121 The figure 
indicates that there is a phase transition at critical = k c (V). For the 
observed volumes k c appears to be small and negative (k c (8000) ~ —0.3 
and k c (32000) ~ —0.1). For k > k c the fractal dimension is close to two in 
accordance with the supposition that similar to the 4 - dimensional case the 
correspondent phase resembles branched polymers. For k < k c X>(8000) ~ 
3.2 and 12(32000) ~ 4. This is in accordance with the expectations that this 
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phase has a singular nature and corresponds to P(oo) = oo. Further we shall 
call these two phases crumpled and elongated respectively. Our results give 
us arguments in favor of view that the nature of these phases is similar to 
that of the phases of the lower dimensional models. We must notice, however, 
that the complete proof has not yet been obtained. 

4. The linear size of the system can be evaluated as Vd. So, we expect 
that in the crumpled phase the diameter d ~ 10 while in the elongated phase 
d ~ 100. These expectations are in accordance with the direct measurements 
of the diameter as well as another parameter called linear extent (see [7j). 
The linear extent £ is defined as the average distance between two simplices 
of the triangulation: 

£ = 7^E<pM> ( 10 ) 

* u,v 

Due to its construction £ should be close to half a diameter. In our mea- 
surements we calculate the linear extent using the slightly different definition 
(that, anyway, should lead to the same result after averaging over the mea- 
surements). Namely, we calculate 

C(s) = - £ < P(s, v) >= i £ R(V R (s) - V R ^(s)). (11) 

v v V R 

Performing the averaging over the measurements we obtain the required value 
of the mean linear extent. Our results on £ are represented in the figure 
Fig. El One can see that the linear size of the manyfold is increasing very 
fast in the elongated phase, while in the crumpled phase it remains almost 
constant (and close to that of 3, 4 and 5 dimensional models j2J). We found 
that the fluctuations of the linear extent almost absent in the crumpled phase 
and are of the order of the mean size of the manyfold in the elongated phase. 



4 Discussion. 

The force between two particles (of masses m) in the classical mechanics is 
equal to 

Gm 2 

F = -% (12) 

where R is the distance between them. Therefore negative gravitation cou- 
pling G and, correspondingly, negative K2 in the four - dimensional model 
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would correspond to repulsion instead of attraction. The same picture is, 
of course, valid for the negative kd_2 in higher - dimensional models. This 
means that the description of gravity in the dynamical triangulation ap- 
proach may appear only for the positive kd-2- The phase transition in the 
pure Euclidean four - dimensional gravity is at K2 ~ 1 and the model exists at 
physical couplings both in the crumpled and in the elongated phases. In the 
ten dimensional model for the observed volumes the phase transition point 
is ~ —0.01. So at physical gravitational couplings the system cannot exist 
in the crumpled phase. 

The fractal dimension at G > is found to be close to 2. This means that 
in order to construct a dynamical triangulation model that corresponds to 
the dynamics of Riemannian manyfold (of Euclidian metrics) of the dimen- 
sion V > 2, it is necessary to add something to the pure Euclidean gravity. 
However, the investigation of the latter seems sensible, because mechanisms 
that appear in this model may play a role in more realistic models. From 
this point of view the most interesting is the mechanism of compactification. 
We started from the model, which has D (^~ 1 ) = 45 local degrees of freedom 
correspondent to the independent components of the metric tensor g. Finally 
we arrive (in the elongated phase) at the system, which has the fractal dimen- 
sion close to two. We expect that similar to the lower dimensional cases this 
system can be effectively approximated by a tree graph. In the correspondent 
branched polymer model |2] the lengths of the linear pieces of the graph and 
the way of their gluing together are the dynamical variables. What does then 
happen with the remaining degrees of freedom? They correspond to the small 
fluctuations of the pieces of the original manyfold represented by the links 
of the graph. Those fluctuations, obviously, become the internal degrees of 
freedom living on the polymer. This dimensional reduction should, in prin- 
cipal, resemble the dimensional reduction of the Kaluza - Klein models [TT] . 
However, the concrete realization of the mechanism has to be investigated. 
We expect that more realistic realizations of the dynamical triangulation ap- 
proach should possess the dimensional reduction, which may be similar to 
that of the pure Euclidean gravity. 

The investigation of the 5 - dimensional model (see j3]) indicates that it 
has a complicated phase structure. Namely, at k,d-2 < —5 several different 
vacua appear. And the system possesses tunnelling between them. In our 
research of the 10 - dimensional gravity we restricted ourselves by the ranges 
of Kg G (—0.1,0.1), where the phase transition between crumpled and elon- 
gated phases is observed. However, taking into account the above mentioned 
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property of the 5D model we expect that the phase structure of the consid- 
ered model is, probably, not limited by the observed two phases. Another 
open question is the order of the observed phase transition. Although we 
have some indications that it is of the first order, the complete investigation 
of the subject has not been performed. 

5 Conclusions. 

In this paper we report our results on the numerical investigation of the 10 - 
dimensional Euclidean quantum gravity in the framework of the dynamical 
triangulation approach. Our summary is as follows. 

1. The considered model contains the phases that we call crumpled and 
elongated. The arguments in favor of view that those phases resemble the 
correspondent phases of lower dimensional models were obtained. 

2. For the observed volumes the phase transition between the mentioned 
phases corresponds to small negative values of k 8 . Therefore at the "physi- 
cal" positive gravitational coupling G ~ the model does not contain the 
crumpled phase contrary to the lower dimensional cases, where critical value 
of kd-2 is positive. This also means that the natural state of the collection 
of 10 D simplices (that appears when the Einstein - Hilbert action is turned 
off and k 8 = 0) belongs to the elongated phase (at least for the observed 
volumes). This is in the contrast with the behavior of the low dimensional 
models: at D = 3,4 and = the systems exist in the crumpled phase. 

3. The average scalar curvature appears to be positive at the observed 
gravitational couplings (k$ e (—0.1,0.1)) in contrast with, say, the 4 - di- 
mensional case, where it becomes negative at k 2 close to 0. 

4. The fractal dimension of the model is close to 2 in the elongated phase 
for both considered volumes (V = 8000 and V = 32000). In the crumpled 
phase the fractal dimension appears to be ~ 3.2 for V = 8000 and V ~ 4 for 
V = 32000. 

5. The linear size of the system in the crumpled phase appears to be 
close to that of the lower dimensional models (the linear extent is ~ 10). In 
the elongated phase the mean linear extent and its fluctuations appear to be 
sufficiently larger than in the crumpled phase. 

6. It is expected that the considered phase transition is of the first order. 
However, the investigation of this subject has not been performed. So the 
question about the order of the phase transition remains open. 
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Figure 1: The curvature. 
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Figure 2: The fractal dimension. Data for V = 8000 is represented by the 
points and the solid line. Data for V = 32000 is represented by the circles. 
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Figure 3: The linear extent. Data for V = 8000 is represented by the points 
and the solid line. Data for V = 32000 is represented by the circles and the 
dashed line. 
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